21 cm radiation - a new probe of variation in the fine structure constant 
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We investigate the effect of variation in the vaiue of the fine structure constant (a) at high redshifts 
(recombination > z > 30 ) on the absorption of the cosmic microwave background (CMB) at 21 cm 
hyperfine transition of the neutral atomic hydrogen. We find that the 21 cm signal is very sensitive 
to the variations in a and it is so far the only probe of the fine structure constant in this redshift 
range. A change in the value of a by 1% changes the mean brightness temperature decrement of 
the CMB due to 21 cm absorption by > 5% over the redshift range z < 50. There is an effect of 
similar magnitude on the amplitude of the fluctuations in the brightness temperature. The redshift 
of maximum absorption also changes by ~ 5%. 



Introduction. — Why do the fundamental constants of 
nature take the values that we measure ? This question 
was posed as a fundamental problem of physics in the last 
century. Dirac [l|, [2| first considered the question of vari- 
ation of these constants with time. The standard model 
does not explain the values of these constants, especially 
the constants determining the strength of the four funda- 
mental forces. These constants can indeed vary naturally, 
though not necessarily, in space as well as time in Grand 
Unified Theories (GUTs) and theories of quantum grav- 
ity 0]. Thus a measurement of variation or a constraint 
on non-variation of the fundamental constants is an im- 
portant probe of new physics beyond the standard model 
and general relativity. Testing the variation of constants 
has become more important in light of current data in- 
dicating the presence of dark energy. Dark energy could 
be a cosmological constant, which fits current data 0], 
but could also be evidence for physics beyond the stan- 
dard model. It is therefore important to find new ways 
to distinguish different models of dark energy. If dark 
energy couples to the standard particle physics it could 
cause variations in the fundamental constants of the stan- 
dard model [j| . Testing for variations in the fundamental 
constants thus also probes the properties of dark energy 
[E 0] and hence is important for any attempt to explain 
dark energy theoretically. 

In this Letter we propose a new probe of variations in 
the fine structure constant on cosmological time scales. 
The most stringent existing constraints on the value of 
the fine structure constant in the early Universe are from 
the measurements of quasar spectra involving fine struc- 
ture transitions. These measurement suggest a change 
| (a< — a)/a| = |<$a/a| < 10~ 5 , where a t is the value of 
a at time t [1, H, EH • For these quasar measurements t 
corresponds to a z < 3.5. If the variation is monotonic 
\5a\ would be larger at higher redshifts. If the variation 
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is not monotonic, then measurements at many different 
redshifts would be required to trace the features in its 
evolution. The constraints at very early times, z > 10, 
are not very stringent and come from the CMB (z ~ 10 3 ) 
and bi g ba ng nucleosyntheis (BBN, z ~ 10 9 — 10 10 ). From 
CMB 0, Qj, \Sa/a\ < 3 - 9% and from BBN El 
|<5a/a| < 6%. There are no constraints for redshifts in 
the range 1000 > z > 10. 

The 21 cm absorption of CMB provides an opportunity 
to constrain a during these "dark ages." Also the 21 cm 
absorption signal is available for a range of redshifts dur- 
ing this period and thus can be a useful probe for trac- 
ing the evolution of a. One advantage of this method 
is that we are measuring the amount of absorption of 
the CMB radiation which depends on the Einstein emis- 
sion/absorption coefficients. As we will see below the 
Einstein coefficients are more sensitive to changes in a 
(Aio oc a 13 ) than the fine structure/hyperfine structure 
splitting itself. 

21 cm radiation from the dark ages. — The 21 cm sig- 
nal from neutral hydrogen after recombination and be- 
fore re-ionization has been investigated by many authors 
EE EE E3, Ell - We refer to 1!) for detailed review. Af- 
ter recombination the radiation temperature goes down 
as 1 + z. The baryons however are prevented from cool- 
ing adiabatically due to the small amount of residual 
electrons which couple the gas to the radiation through 
Thomson scattering. At z ~ 200 this process becomes 
inefficient and the matter decouples thermally from the 
radiation. The hydrogen atoms after recombination are 
in the ground state which is split into a singlet and a 
triplet state due to hyperfine splitting. The occupation 
of the excited triplet state and the lower energy singlet 
state can be described by defining the spin temperature 
by the relation n t /n s = g t j g s e~ T * I Tb , where n t and n B 
are the number densities of atoms in triplet and singlet 
states respectively, gt and g s are the corresponding statis- 
tical weights with gt/g s — 3, T* = 0.068X is the energy 
difference between the two states and equals 21 cm in 
wavelength units and T s is the spin temperature [2pj . 
The evolution equations for the gas temperature T g , 
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radiation temperature T 7 = 2.726(1 + z)K, ionization 
fraction x = n e /nn, where n e is the number density of 
electrons and uh is the total number density of hydrogen 
nuclei, and spin temperature T s can be written as [HI, Hi], 



dT s _ 4T S 2 
HI ~ H (1 + z) \T g 
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where H is the Hubble parameter at redshift z, 
Cio = [Kxo tih + k\ qXiih) , ^uf^ ^ s the collisional de- 
excitation rate from triplet to singlet state for H-H colli- 
sions [IE HE] , K^n is the corresponding cross section for 
e-H collisions (25[, Aiq is the Einstein A coefficient for 
spontaneous transition, ax is the Thomson scattering 
cross section, usb is the Stcphan-Boltzmann constant, 
m e is the mass of electron and c is the speed of light in 
vacuum . xjj e = n-He/n-H, n-He being the number den- 
sity of Helium nuclei. C r , (3 and 012 are defined by the 
equations |1 H, H| 
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where A a = 8irH (1 + z) /A 3 , A Q = 8?r?ic/3Bi and 
£?i = a 2 m e c 2 /2. F = 1.14 is the fudge factor to take 
into account the non-equilibrium among higher energy 
levels of hydrogen [2(| . Helium recombination can be ig- 
nored since it has no effect on the 21 cm signal. Solving 
equations QJ3]) in a given cosmology gives the evolution 
of spin temperature with redshift. The change in the 
brightness temperature of the CMB at the correspond- 
ing wavelength is then given by [13, H3| 
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where u 2 i = k B T+/h ~ U20MHz. The brightness 
temperature is related to the observed intensity by the 
Rayleigh-Jeans formula T\, = I„c 2 /2ksv 2 , where I v is the 
specific intensity and v is the frequency of observation. 

There will also be fluctuations in T s and T& due to in- 
homogeneities in tih and T g which are related to the 
primordial inhomogeneities in the gravitational poten- 
tial (l7l [l8| . The angular power spectrum is then given 
by Ci{z) = {ai m a*i m ), where a\ m are the expansion co- 
efficients in the spherical harmonic expansion of 6T{> = 



Tft — Tb, where Xf, is the mean brightness temperature. 
We follow [TH in our calculations. The baryon power 
spectrum is calculated using CMBFAST [H|. 
Effect of change in a. — A different a during recombi- 
nation affects the CMB power spectrum due to Thom- 
son scattering of photons through equation © and the 
Thomson scattering cross section [2^, 30]. As with CMB 
we ignore the variation in the fudge factor, F, with 
a in Eq. since its effect is negligible. The cru- 

cial point is that for the 21 cm transition, there is the 
following additional dependence on a in equations |[T]). 
(|2|) and 0. The Einstein A Coefficient is given by 
Aw = 647r 4 i/| 1 S , 2i/3ft.c 3 <72, where g 2 is the statistical 



weight of excited state, 5; 



3/?m, /3m = eh/4:Trm e c 



being the Bohr magneton [3l|. Now U21 
and T+ oc ^21, where i?oo is the Rydberg's constant. 
Thus we see that the spontaneous emission coefficient, 
A10 oc v 2 i 0\[ oc a 13 , is a very sensitive function of a. 

The a dependence of the collisional de-excitation rate 
is more complicated. We use ab-initio calculations and 
asymptotic formulae for large separations of potential 
energy curves of the ground state and the first excited 
triplet state of hydrogen molecule to calculate the spin 
change collision cross sections [H, [H, UH, HE HE H3- 
These are nothing but the total energy at a given sep- 
aration in the clamped nuclei approximation or the ex- 
pectation values of the electronic Hamiltonian, H e , which 
has a kinetic energy term (T) and a Coulomb potential 
term (V): 



H e =T + V 



(8) 



A change in the fine structure constant can be treated 
as a perturbation in the Coulomb potential (V oc a). 
Therefore if S is the fractional change in alpha, so that 
ctnew — a{l + 5), then V new = V (1 + 6). Now to first 
order in perturbation theory the expectation value of the 
new Hamiltonian is given by 

(H new ) = (T + V new ) - (H e ) + S(V). (9) 

In above the () denote the expectation value over the 
unperturbed wavefunction everywhere. Also the change 
in the derivative to first order is given by 8dV/dR. We 
estimate dV/dR by fitting a polynomial function to the 
V-R curve from the ab-initio calculations [H, HE HE HE 
HE H3 • Thus from the ab initio potential energy curve at 
unperturbed a, we can construct the first order correc- 
tions when a changes by a small amount. We first calcu- 
late the scattering phase shifts by integrating the partial 
wave equations and use them to calculate the scattering 
cross sections and the rate coefficients (K^ n H ) using the 
standard scattering theory [2E HE HE HE] • 

We checked our code by comparing our cross sections 
for unperturbed a with those calculated by Zygelman 
and Allison & Dalgarno [IE HE] ■ They agree to better 
than 1% for T g > 40K and to 2% for 15K < T g < 4QK, 
The small error that we make at low temperatures is due 
to our ignoring the higher order effects which have been 
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FIG. 1: The fractional change in k^ h with respect to the 
unperturbed value. 
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FIG. 2: Upper panel shows T, with different values of a. 
Bottom panel shows the fractional change (T& (a) — Tf,)/T&. 



taken into account in the calculations in [2J]. We ver- 
ified that this small disagreement at low temperatures 
has a negligible effect on our results by repeating all the 
calculations using the cross sections in [23j, |24( for unper- 
turbed a. Fig. Q] shows k^ h and the fractional change 
in k^q H as a function of temperature for different values 
of a. Changes in k\q due to variation in a has negligible 
effect on Tf, and can be ignored. 

Results and Discussion. — We use ACDM cosmology 
with WMAP3 parameters [i[ in our calculations. Fig. 0] 
shows the observable T, as a function of observed red- 
shift (z bs = V2\ w ' /v bs — 1) where v a \, s is the observed 
frequency today. From Eq. Tf, oc Aiq/v^ oc a 5 giving 
ATb/Tb = 5% for 1% change in a. This is approximately 
the change we see in fig. @] Change in a also changes 
the coupling of T s to T 7 which is opposite to the above 
mentioned a 5 effect. Also z a b s corresponds to different z 
for different values of a causing additional change in the 
signal, ft is clear that the 21 cm brightness temperature 
is a sensitive probe of the variations in the fine structure 
constant. The maximum relative change in the bright- 
ness temperature (ATb/Tb) is > 5% at z Q b s < 50 for a 
1% variation in a. 
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FIG. 3: Upper panel shows the angular power spectrum 
\/l(l + l)Ci/2iv at several rcdshifts. Bottom panel shows 
(y/Ci (5a = -2%) - VC~i)/VWi for the same redshifts. 



Variations in a also affects the power spectrum of 
the spatial fluctuations. This is shown in fig. [5] for 
5 a = —2%. The amplitude of the fluctuations is pro- 
portional to Tb- Thus we expect the amplitude of the 
power spectrum to have similar dependence on a as Tf,. 
There is however also contribution due to fluctuations in 
T s due to inhomogeneities in nu and T g . This additional 
effect causes an increase in the sensitivity to variations in 
a compared to Tf, at z a b s > 100 and a decrease in sensi- 
tivity at z bs < 100. There is also a change in the sign of 
A^/Cl at z bs ~ 40 which is different from Tf,, where this 
occurs at z bs ~ 280. This change in sign is characteristic 
of the a dependence. 

The detectability of the signal is limited by noise due 
to foregrounds. The noise can be expressed in tem- 
perature units for a single aperture telescope as 
Tv = T S ysl 1 e(Avtint) 1 ^ 2 where e is the aperture effi- 
ciency which is close to unity, Av is the bandwidth, tint 
is the integration time and T sys is the system tempera- 
ture which at low frequencies is the temperature of the 
galactic foregrounds ~ 19000K at 22 MHz in a quiet por- 
tion of sky gj. T N ~ 1.4mK for Av ~ AM Hz and 
t; m t = 2000/irs. This means that the sensitivity of a sin- 
gle station of a telescope like LWA Q or LOFAR [H 
can give a constraint on Aa of ~ 0.85%, improving to 
~ 0.3% for the full LWA. The fundamental challenge to 
realizing this measurement is the required precision to 
which foregrounds have to be subtracted. 

Although the foreground removal from the mean sig- 
nal may prove to be impossible, this problem may be 
overcome by using the effect of 5a on the angular power 
spectrum of fluctuations in the 21cm absorption. Sev- 
eral promising foreground removal strate gies have been 
proposed in the literature 0, [4^, Hfl |47| . To ob- 
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tain a similar sub-percent constraint from the angu- 
lar power spectrum requires the higher sensitivity at 
small angular scales (Tjv ~ O.lmK) of a larger tele- 
scope (for example a low frequency equivalent of SKA 
, http://www.skatelescope.org) or much longer observa- 
tion times This is likely to be feasible with future 
technological advancement. 

This new probe is complementary to CMB and BBN 
since it is in a different redshift range and has the poten- 
tial to provide constraints comparable to the CMB ex- 
periment Planck The 21 cm signal is of course also 
sensitive to the cosmological parameters. A 1% change 
in the baryon density, f2&, has a ~ 2% effect on Tb while a 
1% change in the Hubble parameter changes T& by ~ 3%. 
Similar change in the Helium fraction from BBN and 
the matter density, O m , change Tb by < 0.5%. Future 
CMB [48| and large scale structure experiments [45^ [5(| 
will be able to determine these parameters to ~ 1%. As 
seen earlier, the variation in the 21 cm signal due to vari- 
ation in a shows a characteristic dependence on z. This 



is difficult to mimic by changing the cosmological param- 
eters and provide a way to disentangle the two. In the 
power spectrum there will be additional degeneracy due 
to the initial conditions which might be difficult to sep- 
arate. The 21 cm signal is also affected by variations in 
the gravitational constant (G) and the electron to proton 
mass ratio (/x). A complete treatment should consider 
variations in all the constants simultaneously. 

We note that the 21 cm observations can also provide 
a test for the spatial variations of a. This can then be 
used, in principle, to probe the spatial perturbations in 
the dark energy in the framework of theories of quantum 
gravity which relate a to dark energy. 
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APPENDIX 



Plots corresponding to figures 2 and 3 but plotted 
and compared at actual redshift z. They are also avail- 
able as EPAPS Document No. E-PRLTAO-98-004710, 
http : / / www. aip . org/pubservs /epaps . html . 




FIG. 4: Upper panel shows Ti, with different values of a at 
actual redshift z. Bottom panel shows the fractional change 
(T b (a) - T b )/T b . 
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FIG. 5: Upper panel shows the angular power spectrum 
y/l(l + l)Ci/2iv at several redshifts. Bottom panel shows 

(y/Q (5a = -2%) - VC~l)/VCl for the same redshifts, each 
compared at actual redshift z. 



